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Abstract — A general decentralized computational framework 
for set-valued state estimation and prediction for the class of 
systems that accept a hybrid state machine representation is 
considered in this article. The decentralized scheme consists of 
a conjunction of distributed state machines that are specified 
by a decomposition of the external signal space. While this 
is shown to produce, in general, outer approximations of the 
outcomes of the original monolithic state machine, here, specific 
rules for the signal space decomposition are devised by utilizing 
structural properties of the underyling transition relation, 
leading to a recovery of the exact state set results. By applying a 
suitable approximation algorithm, we show that computational 
complexity in the decentralized setting may thereby essentially 
reduce as compared to the centralized estimation scheme. 

I. Introduction 

Set-valued state computation is often used in the analysis 
and synthesis of complex systems. As state sets are thereby 
guaranteed to contain the true state of the system, such 
a computational approach can be efficiently employed for 
the prediction of the system's behavior involving physical 
and measurement uncertainties. For instance, reachability 
analysis for verification of safety specifications is a typical 
application in this context. In hybrid systems, abstraction- 
based approaches naturally lead to a set-valued compu- 
tational framework, see e.g. [1]. Yet, the computational 
complexity remains often prohibitive, which has been an 
impetus for the increasing interest on the decentralized state 
estimation and prediction schemes, particularly in the area of 
discrete event systems with applications to failure detection 
and diagnosis, see e.g. [2], [3]. 

In this paper, a general framework for decentralized 
estimation and prediction for a class of hybrid and dis- 
crete event systems with a finite external signal space is 
considered. Therefore, initially the original signal space is 
decomposed into a finite number of aggregate signal spaces 
of a lower cardinality, which physically may be interpreted 
as substitution of a "high resolution sensor" by a finite set 
of "coarser" ones. Thereby, each of the introduced aggre- 
gate signal spaces invokes a distributed state machine by 
"relabeling" the symbols of the monolithic state machine. 
However, due to the reduction in the measurement resolution, 
state set predictions of individual distributed state machines, 
and as a consequence, of the whole decentralized scheme 
itself, is indeed conservative. To obviate this, the "coarse 
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sensors" must be appropriately designed in the sense that the 
intersection of the individual computation outcomes resulting 
thereof, invariably leads to exact state set estimates. To 
this end, in this article, specific signal space decomposition 
rules are constructed by utilizing the structural properties 
of the transition relation of the underlying monolithic state 
machine. A simple algorithm is developed based on the 
concept of so-called "non-deterministic chains", which rep- 
resent a special class of transition relations featuring inherent 
injective transition functions. Our computational framework 
relates to the work [4] which focuses on distributed set- 
valued state estimation for I/S/O machines. Yet it includes 
additional perspectives of decentralized computation. For 
instance, one of the primers here has been the development 
of the guidelines for the design of the decentralized scheme 
itself. Moreover, our framework covers the larger class of 
state machines that preserve injectivity in the state transi- 
tion function such as non-deterministic automata or I/S/- 
representations with singleton output maps. Recently, this 
approach has been further generalized in [5] to the class of 
systems involving non-injective state transition functions. 

The remainder of the paper is organized as follows. In 
Section|ll] we recall a few basic concepts from the behavioral 
system theory. In addition, a recursive algorithm for the 
computation of the set-valued state estimates and predictions 



is presented. Section III represetns the core of the work. 
Here, we define the decentralized computational scheme, 
and introduce external signal space decomposition using 
specific aggregation functions. In particular, the concept 
of non-deterministic chains and the corresponding main 
algorithm for the signal space decomposition are presented. 



In Section IV the results are employed in a decentralized 
estimation setup using ^-complete approximations. In addi- 
tion, a comparison of the space/time complexity between the 
decentralized and monolithic scheme is shortly discussed. 
Finally, an illustrative numerical example for decentralized 
set-valued state estimation is included. 

We use the following notation: capital letters denote signal 
spaces, e.g. X, U, Y, W represent the state space, the input 
and output spaces, and the external signal space, respectively. 
The corresponding elements are denoted by greek lowercase 
letters, e.g. W ~ {uji, . . . ,uj„i}. We consider the discrete 
time domain, hence signals, which are denoted by lowercase 
letters, are sequences of symbols from the appropriate signal 
space, e.g. w : Nq ^ W represents the external signal. The 
restriction of a signal to an interval [r, t], T,t E No, < 
T < t, is denoted by ■\[T,t]^ ^-8- w|[r,t] = uj{t) . . .w{t). The 
space of the finite sequences (strings) w|[r,t] will be denoted 



as M^I"^'*] — VF*""^"*"^. The string iu|[T-.t] will be considered 
as an element of Vl^I'^ *!, i.e. they will be represented by an 
(t—r+1) -tuple ordered by the time parameter. Finally, let / be 
an arbitrary function / defined on some domain X. If S is a 
subset of X, we will use the convention /(5) :— Uxesfix)- 
For singleton sets we avoid usage of brackets. 

II. Preliminaries 

A. Systems & realizations 

This section provides basic system and realization con- 
cepts from the behavioral perspective, see e.g. [6]. A dy- 
namical system E is defined as a triple (T, W,B), with time 
axis T C M, the external signal space W, and the behavior 
B C W'^, where W'^ ^ {w : T ^ W}. Throughout 
this article, the discussion is confined to discrete-time sys- 
tems, hence T = Nq. The external signal space is finite: 
W = {wi, . . . , Wm}. B then represents a set of sequences 
w : No — ^ which are compatible with the dynamics 
of the system E. A dynamical system E = {Hq,W^B) 
is said to be time invariant if aB C B, where a is the 
backwards time shift operator: aw{t) := w{t + 1), < G Nq, 
and aB := {crw; w G B}; for t > 1, cr^ := ctct^^^ . 

A state machine is a tuple P — (X, W, A, Xq) where X 
denotes the state space, W the external signal space, A C 
X X VF X X the ti-ansition relation, and Xq C X the initial 
state set. If X = M" x D, where n e Nq and e N < oo, 
P is a hybrid state machine; for n = 0, P is a finite state 
machine. Throughout the paper, P is assumed to be non- 
blocking, that is for all ^ e X there exists a. uj ^ W such 
that (^,a;,^') G A. Furthermore, we assume = X. Then, 
aB = B, implying that E is time-invariant. 

For systems exhibiting an input/output structure, the exter- 
nal signal space W can be decomposed as W = U xY, with 
U and Y being the sets of input and output symbols. Then, 
P = {X,U X Y,A,Xq), is said to be an I/S/- machine if 
for each state (, E X and each /i G [/, there exists a ly E Y 
and a ^' e X, such that {£_,{fi,iy),C) E A. If ly and C 
are unique for all ^ G X and fi E U, P is said to be an 
I/S/O machine. Note that I/S/O machines are deterministic 
by definition. A state machine P = {X,W, A, Xq) induces 
a state space system Es = (No, 14^ x X,Bs), where Bs is 
referred to as the full behavior, and is defined as 

Bs:^{iw,x);{x{t),wit),xit+'^))'^^,tENo,xoEXo}. (1) 

The external behavior Sex of Es is then defined to be 
the projection of Bs onto W^", that is Sex := VwBs = 
{w,3x E X^", {w,x) E Bs}. Finally, a state machine P ~ 
{X, W, A, Xq) with induced external behavior Sex is said to 
be a realization of a dynamical system E = {Nq, W,B) if 
Sex = B. This will be denoted by P ^ Y.. 

B. State set estimation & prediction 

Let Ss and Sex be the induced full and external behavior 
of the state machine P — {X, W,A,X), respectively. Define 
state sets compatible to an external string w|[r,ti E l^I"^'*! at 



the time instants t and t + 1 as 

x{w\[T,t]}-={S.]^{w'x)EBs,x{'t)^i^w'\[r^t]^w\[r,t\}, (2a) 

P{w\[T,t\)-H^]^{w'x)EBs,x{t+l)=^,w'\[r.t]^w\[r^]}- (2b) 

Both, X and p, are families of set-valued functions — 
2^ parametrized by a restriction interval [T,t\. Note that 
w\[T,t\ G 'Sexlfr,*] ^ x{w\[T,t]) ^ 0. In general, more 
information about the past leads to more accurate state 
estimates. This fact is reflected by the following proposition. 

Proposition 1: Consider a machine P ~ {X,W, A, X) 
with induced external behavior Sex, and let w E Sex- Then 

x{w\[o,t]) C x{w\[i,t]) C . . . C x(w|[t^f]), (3a) 

p(w|[o,t]) ^ Pi'w\[i,t]) ^ • • • C p{w\[t4])- (3b) 
Proof: Introduce the set of behaviors Br ■— {w' E 
Sex; "i^'IfT,*] — ^\[r,t]}' which contains all sequences in Sex 
that share the same restriction wlj^ t]. Then, by definition 
<|2a]), xiM[T.t]) = Xi^r) ■■= Uu,'gB^x(w'|[r,t])- It is obvious 
that So C Si C . . . C St, which impHes xi^o) Q xi^i) C 



• • • ^ xi^t)- This proves (3ai. Equation (3b i follow by same 
lines of argument. ■ 
We introduce now few one-step prediction expressions, 
which will reveal an iterative computation procedure for 
the state sets in ( |2a| and ( |2b] i. Therefore, introduce the 
parametrized state transition function p^i : X 2^ as 



(4a) 



For fl C W and ^ C X in accordance with the adopted 
notation convention we define 



(4b) 



The predicted states in X, resulting from the occurrence of 
the symbol cj G W, are then computed by 



p{uj) = Pcu(x(^)) 
where according to (j2^ 

xH = U;3e',(e,^,e')e A}, 
with ljj — w|[t,t]- For sequences, we get analogously 

PHlr.t]) = Pw{t)ix{w\[r,t]))- 

Moreover, by definition ( |2a| l 

X{w\[r,t]) = p{w\[rj-l]) n X{w{t)), 



(4c) 



(4d) 



(4e) 



(4f) 

which along with ( [4e] i reveals a recursive structure in com- 
puting x{w\[r,t]) and p{w\[r,t])- 

III. Decentralized computation scheme 

A. Signal space decomposition 

Consider the external signal space W, as defined in 
Section [n] and introduce a finite set of aggregation functions 

Ak-.W^Vk, fcG {l,2,...,p} (5) 



where \Vk\ < \W\. The functions (jsjl are required to fulfill 
the following resolvability or consistency condition 

nLi-4^'(AH)=^ {u^ew) (6) 

where the inverse mapping A].^ : Vfe — > 2*^, is defined as 



(7) 



Due to consistency, each symbol w e Vt^ is uniquely resolved 
by an ordered p-tuple (6'i, . . . , Op), where 9k — y^/c(w), k € 
{1, . . . We refer to this as a decomposition of the original 
signal space W into p signal spaces, and write 



X ^2 X . . . X K. 



(8) 



In general, the opposite does not hold: not every p-tuple in 
Vi y. V2 y- ■ ■ ■ y. Vp will be associated with a symbol in W . 

Now extend the definition of aggregation functions to 
Ak ■ W^'^'*^ — >■ Vj!^'*^, by a symbolwise mapping. That is, 

Vk\lr,t] = A(w|[r,t]) if Vkil) = Akiw{l)) {t < I < t). 

Then, it is clear that the resolvability condition (|6]) carries 
over to strings, as well, that is 

ri^^Ak\Ak{w\^r.,i\)) = M[rA {weW^^>). (9) 
B. Distributed state machines 

Having introduced the signal spaces Vk, k £ {1, . . . 
each is now associated with a distributed state machine Pk — 
(X, Vfe, Afc,X), where AfeCXxV^xX is defined as 



Afe - {{i, Ok. 0;^^^Al\ek), (e, A}. 



(10) 



The original state machine P = (X, W,A,X) is referred 
to as the monolithic state machine. From ( fTO] ) it follows for 
the full and external behavior of a machine Pk- B^.k '■— 
{{vk,x);vk = Ak{w), {w,x) e Ss} andSex,fc = A(Sex) := 
{Ak{w)\w G Bex}, respectively. The definitions for the 
estimation and prediction functions Xk ■ ^i^'*' ~^ 2^ and 
Pk : Vj!'"'*' — ?> 2^, k e {1, . . . ,p}, are analogous to those in 



( 2apb I for the monolithic machine P. Equivalently stated: 

Pk = poA^\ Xk^XoA^\ (11) 

For instance, 

PkiVk\[rJ]) = P{Ak'^{Vk\[r,t])) = U,= lP(w^fci|[r,t]), (12) 

for some Sk E N. The external behavior B^x.k of the 
machine Pk is a coarse approximation of the "monolithic" 
behavior Sex, in that w\[^^t] e A^^{vk\[r,t]) if Vk = Ak{w), 
for any w e Sex- However, the resolvability condition 
r^l-^iA'i^^ {vk\[rA) ~ M[T.t] suggests a computation scheme 
with p distributed state machines Pk including an intersection 
of the respective outcomes. Therefore, consider a string 
w\[T,t] corresponding to the p-tuple {vi\[r,t], ■ ■ ■ ,Vp\[r,t])- 
Then, in general, it follows 

f^l=lXkivk\[r,t]) = (^l^lXiAk^ivk\[r,t])) 
^X(nLl-4fc'KI[r,t])) 

= x(w|[r,ti), (13a) 



where we use the fact: 0(Afi) n ^(Afa) 3 (j){Mi n M2), and 
the consistency condition. Similarly, 



:lPk{vk\[T,t]) 2 Piw\[r,t])- 



(13b) 



As a consequence, the parallel computation scheme pro- 
vides, in general, an overapproximation of the outcomes of 
the original monolithic state machine. However, for certain 
classes of transition relations A, specific consistent functions 
Ak, k G {l,...,p} can be constructed, which lead to 
exact computation results in the decentralized scheme. The 
basic concept which we use in constructing such aggregation 
functions is that of "non-deterministic chains". 

C. Non-deterministic chains 

Definition 1: Consider a state machine P~{X, W, A, X), 
and let il C W. A transition subrelation (5 C A is said to 
be a non-deterministic chain over il if for all ^ G X and 

(i) (e, Lo',(') e S, (e, Lj", e") eS^u' = uj", 
ill) (e',c.',C)eJ,(e",^",Oe^^e' = r- 

The transition relation 5 can naturally be associated with the 
functions x'^ ■ ^ ^ "2-^ , p'^ ■ — 2^ , and the transition 
function p"^ : x(^) ~^ 2^, defined as 



X 



xUh p 



p\ 



n, P 



Pn- 



(14) 



Then, (i) can be equivalently restated as x'^i'^') ^ x'^i^") 
0^0;' = w", while (ii) is equivalent to p^iC) n p^iC) ^ 
= leading to the following statement. 

Proposition 2: A transition relation (5 C A is a non- 
deterministic chain if and only if and p'^ are absolutely 
injective set-valued maps. 

In the sequel, we introduce a systematic method for signal 
space decomposition which results from partitioning the 
transition relation A into a finite set of non-deterministic 
chains. To this end, for a given monolithic machine P = 
{X, W,A,X), suppose that partitionings 



W 



and A 



(15) 



exist, such that each transition subrelation Sj C A over flj 
represents a non-deterministic chain. Then, we claim that 
the state machine P is chain-decomposable . The following 
example illustrates this idea. 

Example 1: Consider an I/S/- machine P = {X, U x 
Y,A,X) with singleton output maps. Let U = {pj;j = 
l,...,r}, and introduce the partitioning W — U xY = 
Uj=i%, where ilj := {pj} x Y. This induces a partition- 
ing of the transition relation A = U^^iSj. By definition, 
functions f:XxU^2^ and h:XxU^Y exist, such that 

'^=H^,Pi)- (16) 

Then, each state £, E X can be associated with a unique sym- 
bol pair {pj,i') E flj. Hence, {i) in Definition fl] is fulfilled. 
Define further p^j : X 2^ as /5^(C) := f{pj!t), and let it 
be absolutely injective. Then, Sj is a non-deterministic chain 
for all j E {1,. .. ,r}. ■ 



Next, introduce a consistent signal space decomposition 
nTAl for flj, j e {l,...,r} 



Corollary 1 (I/S/-): Consider the class of I/S/- realizations 



as discussed in Section 



X V, 



3,P- 



(17) 



Notice that such decompositions invariably yield a consistent 
decomposition for the original signal space W ~^ViX ■ ■ ■ x 
Vp if, e.g. 



(18) 



The state machines Pk, k E {!,..., p} are then easily 



constructed using the procedure described in Section III 



IbI We now want to show that the intersection of their 
estimates and predictions produces the exact outcomes of 
the underlying monolithic state machine P. 

Referring to ( [T4] l and Proposition [2] it is important to keep 
in mind, that the absolute injectivity of 



(19) 



is, per construction, preserved in all disjoint subspaces ilj. 
Note also that for any t e Nq, all elements of the inverse 
mapping of a symbol 9k E Vk belong to the same subspace 
rij for each k E {!,..., p} and some j E {l,...,r}. 
Now, fix a string w|[r.t] and consider the corresponding tuple 



[T- t], . . . ,Wp|[T.t]) in the decentralized scheme. Then, 

f^l=lPk{vk\[r,t]) = n^=i ^''l!:=l PiWklJ[r,t]) 



U 



U 



=1 ' 'fc: 



-.iPi^kl 



wit)) 



where Sk < Sk is the number of all strings in VFi^ 'i of the 
form WkiJ[r,t~i]w{t), k E {1, • ■ ■ ,p}- In the 3'"'' Une we took 
advantage of the consistency condition and the fact that pj : 
ilj — > 2^ is absolutely injective. All the strings A^^{vk\[T.t]) 
that do not end with w{t) are then neglectable. Moreover, 
using the recursive formula ( |4^ , and the injectivity of the 
transition function p^, 

^l=lPk{vk\lr.t]) = 



= u 



h=l 



^1=1 nLiP,' [x'M^)) n p{wki,\[r.t-i])] 
x'^iwit)) n (uf^^^i • • • uf;^i nLiPKiJ[r,t-i]) 
x-(«^W)n(nLiPfcKI[r,t-i]))]. 



As a consequence, a recursive expression is obtained, which 
transfers the computation task from the interval [r, t] to [r, t— 
1]. Repeating the recursion until [r, r], leads to an expression 
on the right-hand side equal to ^j), which is the proof 

of the following main statement. 

Theorem 1: Let P={X, W, A, X) be chain decomposable. 
Then, a decentralized scheme including the state machines 
Pk, k E {!,..., p} induced by ( jTs] ) provides exact state 
estimates and predictions, that is 

<^l=iXk{vk\[T,t]) ^xMlT,t]), (20a) 

<^k=lPk{'>^k\lT,t]) = p{w\[r,t])- (20b) 

Following the discussion in Example 1 for I/S/- state ma- 
chines, we further conclude. 



P ^ {X,U y. F,A,X), which fulfills (16l. Introduce a 



partitioning of A as described in Example 1. Then, a de- 
centralized scheme built upon any consistent decomposition 
with Vk — U y. Ak{Y), k E {l,...,p}, provides exact 
computation results of the form ( |20a| l and ( |20b[ ). 

Example 2: In order to illustrate the proposed method, con- 
sider the automaton in Fig.[T] Its external signal space can be 
partitioned W — Q,i U VI2 with 
VLi = {ai, 61, ci, rfi} and VL2 = 
{a2, 62, C2, ^2}- As indicated by the 
solid and dashed lines, the cor- 
responding transition relations 5i 
and S2 are both non-deterministic 
chains. According to the previous 
elaborations, any consistent set of 
aggregation functions can be ap- 
plied on rii and il2. E.g. a particular 
signal space decomposition results 
from 




Fig. 1 .Finite state machine. 



Via = with el ^ {ai,h}, d\ ^ {ci, di}, 

^1,2 = with B\ ^ {ai,ci},0i ^ {61, di}, 

= Q\} with ()\ ^ {aM, Q\ ^ {C2, d2}, 
^^2,2 = {QlA} with B\ ^ {a2,C2},e| ^ {hiAi}- 

The corresponding decomposition reads W V\ x V2, 
where V\ = V\_\ U V2,i and V2 = Vi^2 U V2,2> leading to 
the distributed machines Pi — {X,Vi, Ai, X) and P2 = 
{X, V2, A2, X). Now assume that the original system accepts 
a string, saywlf^- t] — 0162- The corresponding estimate of 
the monolithic machine is x(ai62) = ^2- The distributed 
machines measure accordingly the strings Wi|[o,i] = ^16*1 
and W2|[o.i] = ^^2^2' providing the estimates xi(^i^i) = C2 
and X2(^2^2) ^ {^2,^3}, respectively. The decentralized 
estimate is thus given by x(^i^i) ^ x(^2^2) — "^2, which is 
exactly the same outcome obtained by the monolithic state 
machine P = {X,W,A,X). According to Theorem[T] this 
must hold for all strings accepted by P. ■ 

IV. Decentralized estimation using 

i!-COMPLETE APPROXIMATION 

This section addresses the set-valued state estimation for 
time-invariant systems S = (Nq, W,B) using the ^-complete 
approximation algorithm. Note that E is said to be £- 
complete, [7], if 

w E B <^ a*w\it^t+i] l3\[oj], (<eNo). (21) 

A. i-Complete approximation 

Consider a time-invariant system E. The model = 
{No,W, Bi), ^ e N, is a strongest ^-complete approximation 
of S if (i) Bi is ^-complete; (ii) Bi 2 B, and B'^ 2 B, 
B'l being f-complete ^ B'^ 3 B^. We consider here the 
realization algorithm for the strongest ^-complete abstraction 
Ef as follows (see also [4]). 

Definition 2: The state machine P( ~ {Ze,W, Ae, Zq) is 
a realization of E^ with 



(a) Zo - W; 

(b) Zt ■= uf^lW^^ where W = w^^w^^ . . .u;""; 

(c) transition uf^p^f ^ x x Z^, defined by: 

(w I [o,r-i], w(r) , u;| [o,r] ) : w^l [o,r] e I [o,r] , 1 < < ^} , 

^£ :={ I [o.£-i] , , w I ) : w I [o,^] S | [o/] }• 
The ^-complete representation keeps track of the system's 
past trajectories with a sUding time window of length I. The 
resulting behavior Bi includes the sequences compatible with 
E, while So D 61 D • • • D e holds. Using Definition [S] the 
state C G ^£ of P( at a time instant t reads 



l[o,t] 



if < t < 



(22) 



w|[t_^?+i,t] if t>£. 

For t < £, the state <^ is determined by the whole signal 
string, whereas for t > £ only by its suffix of length £. 
Therefore, it is evident that each state C,{t) can be associated 
with unique estimates xiCit)) ^ defined as 



xm) 



(23) 



X(^l[0,t]) ifO<<<^, 
X{w\lt~e+i,t]) if t>£. 

As a consequence of Proposition [T] in the latter case, due 
to x{w\[o,t]) C x{M[t-e-i,t]), the state sets attached to C, 
provide a coarse - yet instantaneous - estimation of the 
system realization P = {X, W, A, X) compatible with the 
measurement w|[o,t] for t > £. Another consequence of 
Proposition [T] suggests that the estimation accuracy can be 
improved by increasing £. Indeed, given two approximation 
automata Pg-^ and Pg^ with £2 > £1, the estimates for 
t > £2 correlate as x{w\[t~t2+i.t]) ^ x{w\[t~ii+i.t])- 
Clearly, the estimation accuracy is improved, as Pi^ stores 
a larger content of information {i.e. a longer suffix) on 
wljo^t]. A major drawback is, however, that the number of 
states and transitions of the ^-complete automaton increases 
exponentially with increasing £, which provides a motivation 
for the decentralized approach. 

For the distributed state machines P^. = (X, V^, Afc,X), 
A: = 1, . . . ,p, as introduced in Section [Ill-B[ the ^-complete 
automata Pk,i = (^fc/, V/c, ^k.i, -^o.fc) are obtained from the 
algorithm in Definition |2] Then, Zk^i includes states 

/ Wfe|[o.t] if Q<t<£, 
Vk\[t-i,t\ if t>£, 



Ck{t) 



(24) 



which store the estimates 



XkiCkit)) :-- 



Xkivk\io,t]) 

Xkivk\lt-Lt]] 



if 0<t<£, 
if t> £. 



(25) 



For a given string w|[7-.t] = (fi| 



\lr,t],---,Vp\[T,t] 



), the mono- 



lithic approximation P( and each Pkj assume unique states 
C ^ Zg and Cfc e Zk^e, as defined above. Then, in accordance 
with the elaboration in Section [nil 

nLiXfc(C.W) 2 xm) or nLi XkiCkit)) = x(CW)- 

B. Space/time complexity 

Following the discussion in the last section, it is obvi- 
ous that the complexity of the monolithic and distributed 
estimators is directly related to the number of states in the 



corresponding implementation with ^-complete automata. In 
accordance with ( 22 1 and ( 24 1, the number \Zk e\ of the states 



of automata with a memory depth £ refers to the number of 
substrings of length £. For instance, for a distributed automa- 
ton Pkj = iZk,i,Vk,Ak,Zo,ky- \ZkA < ELol^kt For 
comparison purposes, here, the maximal number of possible 
states in the automata implementation is considered. That 
is, the numbers X]fc=i Y^i=o l^fcP (decentralized implemen- 
tation) and Yfi=o \ (monolithic implementation) need to 
be compared. As, by definition, \ Vk\ < it is obvious that 
for a sufficiently large £, the decentralized setup requires less 
memory space than the monolithic implementation. 

The memory requirements depend also significantly on 
the amount of the stored data Xfc(C)- For instance, for 
a finite state machine, a finite number of sets has to be 
stored. Similarly, for switched linear systems, the estimation 
sets are polytopes, which, again, can be stored as a finite 
number of sets representing the poly tope vertices. Due to the 
exponential growth of the number of states with increasing 
£, the memory requirements can be significantly reduced by 
the decentralized approach at the price of additional online 
effort for the computation of set intersections. For certain 
classes of hybrid systems this operation can be implemented 
efficiently (such as polytope intersection in the case of 
switched linear systems). For other classes, the estimation 
sets can be offline overapproximated, e.g. by polyhedral 
methods, in order to obtain a computationally efficient online 
intersection operation. In general, a significant reduction in 
the offline computational time and required memory space 
is expected, as exemplified in the following section. 




Fig. 2. Two-tank system. 

C. Numerical example 

The plant consists of two tanks with corresponding water 
levels xi and X2, inflows ui and U2, and outflow parameters 
fli and 02, that are connected by a pipe with a flow constant 
b (Fig.|2]i. The dynamics can be described by the difference 
equations: 



x{t+l) = 



1—ai — b b 
b 1— & 



x{t) 



1 
1 



u{t), (26) 



where t e Nq, x ^ [xi,X2]'^, u = [^1,1/2]^, y = [^1,2/2]^, 
and y{t) = x{t). In the following, the simulations will be 
performed for a set of input symbols C/j ={/ij,/i|,/i^}, where 
Idj, plj represent the inputs uj = 1, uj — 7, Uj — 14, 
respectively, with j E {1,2}. Analogously, Yj = {v}j ,v^}, 
where v^, i/j, encode the measurements yj e [0, 10), yj G 
[10, 20), yj £ [20, 30], respectively. Furthermore oi = 02 = 
0.35 and &=0.25. 

Note that this system assumes an infinite-dimensional 
I/S/O state machine representation P — (X, M^, A,X), 



where the external signal space is given hy W = Ui x U2 ^ 
Yi X ¥2- Hence, in accordance with Example 1, only the 
output space is to be aggregated. For a decentralized setting 
including two state-machines, the resulting signal spaces of 
Pi^e and P2^e are given by Vi ^ Ui x U2 x Ai{Yi x Y2) and 
V2 = U1XU2X A2{Yi X Y2), respectively. Thereby, we adopt 
the mappings Ai and A2 according to {(^,%)} 1— > 9\ and 
^^ 6l|, i e {1,2,3}. As a result, |Fi| + |V^>| =54, 
while \W\ — 81, indicating a reduction in the computational 
complexity (see below). 

The i?-complete approximations of monolithic and dis- 
tributed estimators with ^ = 2 are simulated for < ^ < 
7 with the input sequence u{t) as depicted in Fig.|3] and 
x(0) = [0,0]^. Obviously, the reachable sets for an arbitrary 
sequence of input and measurement symbol pairs is repre- 
sented by polytopes with a finite number of edges. As system 
(26 1 fulfills the conditions of Corollary [T[ the decentralized 
estimation scheme must provide the exact outcome of the 
monolithic estimator. For instance, following ( |22] l and (|24]l, 
at the time instant t ~ 2, the monolithic and distributed 
estimators reach the states 

C(2) = {{fil,fil,iyl,iy2), {fil,fj.l,i^l,i'2), 
Ci(2) = ii^il^ilel), (plfil el), ifilf^l el)), 
C2i2)^{i^l^ilel),ifil,plel),{plt,li 



91)), 



respectively, and the intersection of xi(Ci(2)) (the polygon 
with dashed edges) and X2(C2(2)) (the polygon with dotted 
edges), is indeed equal to the set provided by the monolithic 
estimator (the gray area indicated by x(C(2))), see Fig.|3] 
Note that the true state of the system denoted by small filled 
circles is always contained in the corresponding state set. 

Table |l] indicates that complexity of the estimators is 
strongly related to the number of states \Zi \ of the ^-complete 
automata used for realization of the approximations and the 
overall number of vertices of the associated sets x(C) given 
by n^. Observe that, compared to the monolithic estimator, 
akeady for £ — 2, the decentralized setup with two estimators 
requires a lower offline computational cost \Ze\, and amount 
of memory n^. 



I 






Pi,e & P2,e 




\Zi>\ 


\Zi\ 


T,\Zt\ Sn^ 


2 


172 993 


54 307 


108 617 


3 


2260 17743 


703 5509 


1406 11065 



TABLE 1 
Complexity analysis. 

V. Conclusions 

A general decentralized framework for set-valued state 
estimation and prediction for hybrid state machines has 
been discussed in this article. The outcome of the decen- 
tralized scheme is computed as the intersection of the sets 
provided by the individual distributed state machines. The 
distributed state machines are defined as abstractions of the 
monolithic one, as specified by an appropriate decomposition 
of the original external signal space. In general, we show 
that decentralized schemes provide outer approximates of 



n(6) 
u{7) = [1 7] 

X2(C2(2))^^ 

/ 


= [171^ - ■■ ■■ 

^ / //^=VA 


1 j—^ 

,^f0) = [7 7]^ 


V \i(Ci(2)) 
[7 7]\| 

Xi(Ci(2))nx2{C2(2)) = x(C(2)) 



20 

Fig. 3. Set-valued state estimation using ^-complete approximation (£ = 2). 

the estimates and predictions obtained by the monolithic 
computation. Based on the concept of non-deterministic 
chains, we develop a simple decomposition algorithm for the 
external signal space, invariably leading to exact set-valued 
decentralized state estimation and prediction. Due to the 
smaller cardinality of the constructed external signal spaces 
corresponding to the distributed state machines, a significant 
reduction in the overall space/time computational complexity 
may result. This is illustrated here by applying ^-complete 
approximation. Moreover, advantages in terms of robustness 
and reliability are gained. Indeed, in order to prevail the 
failure of a single processor, the decentralized scheme could 
be extended by a redundant state machine and a failure 
detection for switching off the malfunctioning processor. The 
proposed decentralized framework can be applied in different 
contexts, including data fusion, fault tolerant control, etc. 
Optimal signal space decomposition leading to a minimal 
computational complexity is a problem of interest for the 
future work. A further development of the algorithm based 
on non-deterministic chains can be found in [5]. 
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